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Abstract 
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nuclear Hamiltonian. Up to N'^LO, the Ml operators contain two low-energy constants which 
appear as the coefficients of the non-derivative two-nucleon contact terms. After determining these 
two constants using the experimental values of the magnetic moments of the triton and '^He, we 
have made a totally parameter-free calculation of the hen cross section. The results are in good 
agreement with the data. 
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I. INTRODUCTION 



The radiative capture of a slow neutron on ^He ('^He + n — ^He + 7) , or the hen process, 
is an example of rare situations where the contributions of the single- nucleon (IB) currents 
are strongly suppressed due to the pseudo-orthogonality of the wave functions; viz., since 
the major components of the initial and final state wave functions belong to different repre- 
sentations of symmetry group, they cannot be connected through the r-independent leading 
IB operators. In fact, this suppression is so strong that the meson-exchange-current (MEC) 
"corrections" become comparable to the IB contributions, and furthermore, since the MEC 
and IB terms in this case come with opposite signs, a strong cancellation between them 
occurs, leading to a further drastic suppression of the cross section. The hen process near 
threshold is governed by the Ml operators, and the MEC contributions to them consist of the 
well-known one-pion-exchange part and the short-range part. It is to be noted that the lat- 
ter is not constrained by the symmetries of QCD, and that because of the above-mentioned 
suppression and the cancellation mechanisms, the short-range contributions are crucial even 
for a rough estimation of the hen cross section. Furthermore, the strong suppression of the 
IB matrix elements implies that their values are highly sensitive to the details of the wave 
functions. Therefore, for a precise estimation of the hen cross section, it is imperative to 
have: (i) a reliable method for deriving the relevant MEC operators with a good control 
of short-range physics, and (ii) the accurate wave functions for the initial and final wave 
functions. These requirements make the quantitative estimation of the hep cross section 
highly non-trivial. 

The strong enhancement of the relative importance of MEC in the hen process makes it 
a valuable testing ground for our understanding of the roles of MEC in light nuclei. Apart 
from this point, which is important on its own right, a careful study of the hen process is also 
of great importance in connection with the so-called hep processes, '^He-l-p ^He + z/e + e"*", 
because hep shares all the aforementioned features of hen: the drastic suppression of the IB 
contributions, strong cancellation between the IB and 2B contributions, and the sensitivity 
of the transition amplitude to the details of the nuclear wave functions. The hep is one 
of the proton-burning reactions that take place in the interior of the Sun and, since hep 
produces the highest-energy solar neutrinos among the pp-chain reactions, it has attracted 
much attention in the study of the solar neutrinos; see for a recent review. In a 
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recent paper. Park et al B developed arr effective field theory (EFT) approach, which has 

come to be known as "more- effective EFT" (MEEFT for short) [43, and, with the use of 
MEEFT, they calculated the hep S'-factor with an estimated accuracy of about 15 %. In 
view of the fact that the previous theoretical estimations of the hep S-factor ranged over 
two orders of magnitude [3i], this is a remarkable feat. A direct test of this theoretical result, 
however, is not possible because of the extreme difficulty of measuring the hep cross section. 
Meanwhile, the threshold hen cross section has been measured with a reasonable accuracy: 

n 

(^exp = (54 ± 6) [5] and a^xp = (55 ± 3) fib [6]. Given the close similarity between hen 
and hep, it is expected that comparison between theory and experiment for the hen case 
will offer valuable information on the validity of the theoretical framework employed for the 
hep calculation in This gives an additional motivation for carrying out a detailed study 
of hen. 

Although the theoretical investigation of hen has a long history, the hen cross section 
has never been explained in a satisfactory manner. Towner and Khanna [3] and Wervel- 
man et al. Gj performed shell- model calculations and obtained cr = (14 ~ 125) fib and 
a = (47 ± 18) ab, respectively. Much more elaborate calculations have been performed by 
Carlson et al. |8| and by Schiavilla et al. [9], arriving at a =112 fib and a = 86 fib, respec- 
tively. These works are based on the conventional approach, the so-called standard nuclear 
physics approach (SNPA for short), which consists in the use of phenomenological transition 
operators and phenomenological wave functions. SNPA has been enormously successful in 
correlating and explaining a vast range of electroweak nuclear transitions in light nuclei but, 
from a formal point of view, it has an insufficiency that it is not equipped with a systematic 



way of deriving MEC operators. T 
the approach of Chemtob and Rho 



le MEC operators in SNPA are constructed based on 



MEC operators 

is guided by chiral symmetry and the Ward identities related to the adopted Hamiltonians, 
a certain degree of model dependence is unavoidable. 

In this paper we report on a parameter-free MEEFT calculation for the hen cross section 
at threshold which is based on essentially the same method as adopted in [3].#^ In MEEFT 
the transition operators are derived in a systematic expansion based on HBChPT, and the 



There has been an attempt by Y.H. Song et al. ll| to carry out an MEEFT calculation of hen, but an 



approximate treatment of the scattering wave function in ll| limits its reliability. 
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nuclear matrix elements are obtained by sandwiching these operators between wave functions 
generated from a high-precision phenomenological nuclear Hamiltonian. Mismatch in the 
short-range is overcome by the renormalization procedure for the local operators, according 
to the premise of low-energy EFTs (see below). It thus takes advantage of the systematic 
nature of EFTs and the availability of state-of-art wave functions. The mentioned "formal" 
mismatch may be regarded as a weak point, and also the accurate reproduction of the 
relevant effective-range parameters (ERPs) is not automatically guaranteed in MEEFT; we 
shall come back to these points later in the text. 

In the present work we derive the Ml operators within HBChPT up to N^LO. These 
Ml operators turn out to contain two low-energy constants (LECs), denoted by g^s and g^y, 
which are the coefficients of two-nucleon contact-term operators. These LECs can be fixed 
by requiring that the experimental values of the magnetic moments of the triton and ^He be 
reproduced; this is the same strategy as adopted in Refs. Q, Q, where the Ml properties 
of the A=2 and A=3 systems were studied in MEEFT. A remark is in order here on how 
the short-range contributions are taken into account in MEEFT. The basic premise of EFT 
is that physics pertaining to ranges shorter than the experimentally relevant scale is to be 
lodged in local operators. It means that, provided that a proper renormalization procedure 
is implemented to the coefficients of local operators {i.e., LECs), the renormalization invari- 
ance ensures that the net physical amplitudes be independent of the details of short-range 
physics. We implement the renormalization condition here by adjusting the values of LECs 
{g^s and g^^ for hen and for hep) so as to reproduce a set of known experimental data [ 
/i(^H) and /i(^He) for hen and the tritium-beta-decay rate for hep ]. This matching proce- 
dure should be done for each cutoff value and for each Hamiltonian adopted.*^ Differences 
in short-range contributions calculated for each case shift the values of LECs (which are not 
physical observables) , but the physical amplitudes should remain unaffected if renormaliza- 
tion invariance is to hold. The validity of this scheme can be checked by monitoring the 
stability of the relevant physical observables with respect to changes in the cutoff parameter 
A. It turns out (see below) that, in the present hen case, the inclusion of the local-operator 
(or contact term) contributions reduces the A-dependence by a factor of ~5, demonstrating 
the validity of the adopted renormalization procedure. The residual A-dependence may be 

'^'^ The details of the cutoff parameter A and nuclear Hamiltonians used will be given later in the text. 
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ascribed to higher order contributions. 

On the other hand, for A-hodj systems with A> 3, even a reahstic nuclear Hamiltonian 
often fails to reproduce accurately the ERPs that govern the long-range part of the transition 
matrix elements. If such a mismatch in the long-range region happens, it cannot be cured by 
local operators, which means that one in general cannot expect to have a reliable theoretical 
prediction. As an exception to this general statement, however, we should mention that, 
if a strong correlation between the ERPs and the transition amplitude under consideration 
can be established with high precision, this correlation can be used to drastically reduce the 
model dependence of the calculated transition amplitude [13|. We shall take advantage of 
this feature later in the text (Subsection III A) to obtain as the best estimate for the threshold 
hen cross section: a = (49.4 ± 8.5) fib (for the AV18+UIX potential), and (44.4 ± 6.7) fib 
(for the I-N3L0+UIX potential); see eq. (|T7|) . Good agreement of these estimates with the 
experimental value of the hen cross section gives sound support to the previous MEEFT 
calculation for the hep S'-factor [3f\. 

We wish to emphasize that the present work is the first calculation of hen that employs 
nuclear wave functions which are fully consistent with a given realistic nuclear Hamiltonian; 
these wave functions are numerically exact solutions to the Faddeev-Yakubovski equations 
in configuration space for a specified realistic nuclear Hamiltonian. It is to be noted that 



Big, 



ll| contained some artificial 



hitherto even the most advanced realistic calculations 
adjustments to suppress the ^He(n,p)^H channel in the scattering wave functions. 

This paper is organized as follows. In the next section we explain the formalisms used 
to derive the Ml operators and to obtain the four-body nuclear wave functions. Section 3 
gives the numerical results and analysis. In the final section the implication of our work is 
discussed. 



II. FORMALISM 

A. Electromagnetic current and capture cross section 

We derive here the electromagnetic (EM) currents and Ml from HBChPT up to N^LO in 
chiral order counting. HBChPT contains nucleons and pions as pertinent degrees of freedom, 
with all the other massive fields integrated out. In HBChPT, the EM currents and Ml 
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operators are expanded systematically with increasing powers of Q/A^, where Q stands for 
a typical momentum scale of the process and/or the pion mass, and ~ 4:7^ f-w ~ m ~ 1 GeV 
is the chiral scale, — 92.4 MeV is the pion decay constant, and m is the nucleon mass. 
We remark that, while the nucleon momentum pi is of the order of Q, its energy (~ p1/m) 
is of the order of Q"^ /m, and consequently the four-momentum of the emitted photon g'^ 
should also be counted as O^Q"^ /m). 

The Ml operators up to N'^LO consist of IB and 2B contributions, while three-body 
operators in our counting scheme enter only at N^LO or higher orders. In the center-of- 
mass frame, each current can be written in the form of 

Jem = e-^'-^iiq X j^+jc), (1) 

where q is the momentum carried out by the photon, X = Vi for the IB current pertaining 
to the 2-th nucleon, while X = (r^ + rk)/2 for the 2B currents involving the nucleons j and 
k. The resulting Ml operator in multipole expansion reads 



(2) 



M 



where the subscript M is the magnetic quantum number, g = |q| = 20.578 MeV and jn{x) 
is the spherical Bessel function; the sum runs over all nucleons for IB current and over all 
pairs of nucleons for the 2B current. 

In terms of the hen cross section at thermal energy is given by 

c / q V f Q \ , ,,,2 



with 

A<^(«4t"^»|M„(<,)|9'ii;f-'). (4) 

where a is the fine structure constant, and Vrei = 2200 m/s is the velocity of a thermal 
neutron. 



The explicit forms of the Ml operators up to N^LO are given in our recent papers(l2l. Il3|. 



so we present here only a brief summary of them. The IB Ml operators including the 
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relativistic corrections read 



2rn? / 2m^ 



p1 \ w{2ni - Qi) 



Qifi X « • 1 ^ iVi X (pi X CTj) 



+ -^J2(gri) /ii {3fi Vi- (T,-(Ti)^ I (5) 

where = g is the photon energy, j„(x) = and /ij are the charge and 

magnetic moments of the z-th nucleon, respectively, and is the mean momentum operator 
of the i-th nucleon. 

The two-body Ml operators up to N^LO can be divided into the so/i-one-pion-exchange 
(Itt) term, the vertex corrections to the one-pion exchange (IvrC) term, the two-pion- 
exchange (27r) term, and the contact-term contribution (CT), 

/X2B(g) = [^'7 + (^'/^ + ^'7 + ^S"")] = NLO + N^LO. (6) 

i<j 

In this work we neglect the so-called fixed-current contribution, which is proved to be nu- 
merically negligible In performing Fourier transformation of the 2B operators into 
coordinate space, we introduce a Gaussian regulator with a cutoff A. 

It is to be noted that the contact-term p'?^ contains two low-energy constants (LECs), 
Qis and qav, 



Since the values of LECs are not determined by the symmetry arguments, they need to be 
fixed either by solving QCD at low-energy or by a set of experimental observables that are 
sensitive to these LECs. Since the former is not feasible at present, we resort to the latter. 
Specifically, the values of g^s and g^^v are fixed so as to reproduce the experimental values 
of the triton and '^He magnetic dipole moments, for each nuclear interaction model and for 
each cutoff value. 

B. Faddeev-Yakubovski equations 

Four-nucleon wave functions are obtained by solving the Faddeev-Yakubovski (FY) equa- 

n n 

tions in configuration space [141 . |15| . The FY formalism offers a mathematically rigorous 



description for both continuum and bound states. In this formahsm wave functions are nat- 
urally decomposed into so-called FY amplitudes (FYAs). For A = A systems, there appear 
two types of FYAs, which we refer to as components K and H. The asymptotes of the 
components K and H incorporate 3+1 and 2+2 particle channels, respectively (see FigH]). 
The FYAs are evaluated in isospin formalism, i.e., protons and neutrons are regarded as 



IT'* 

'^12,3 



(j3jz)J®(T3,gT 



34 
2 



FIG. 1: The FY components, K and H. As z — > oo, the K components describe 3+1 particle 
channels, while the H components contain asymptotic states of 2+2 channels, figure [a]. Figure [b] 
shows the j-j coupling scheme used in expanding K and H into partial wave bases. 

the degenerate states with the same mass, which is fixed to fi} jm = 41.47 MeV-fm^. With 
the three-body force Vus = I^i| + V23 + V^^, the FY equations for K = Kf23 and H = iffi 



read 
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E-Ho-Vu-J2 ^S] K = VuiP^ + P') [(1 + Q)K + H] + ^ {V^l + ^3^) v]/, 

i<j J 

J^^Ah = VuP[{l + Q)K + H], 



E-Hn-V, 



12 



(8) 



i<j 



where Vij and are, respectively, the short-ranged part and the Coulomb-dominated long- 
range part of the interaction between the i-th and j-th nucleons. P+ = (P^)^^ = P23P12, 
Q = — P34 and P = P13P24 = P24P13, where is the particle permutation operator. In 
terms of the FYAs, the total wave function of an A = 4 system is given by 



[1 + (1 + P+ + P-)Q] (1 + P+ + P-)K + (1 + P+ + P-)(1 + P)H. 



(9) 



We expand K and H in terms of the tripolar harmonics Y^'lxi, yi, Zi), which comprise spins 
and isospins of the nucleons as well as angular momentum variables, 



-Y^ (Xj, Iji, Zi) , 



(10) 



where $ stands for either K or H, and the subscript i denotes the particle-grouping class 
(among the four nucleons). We note that the total angular momentum and its projection, 
parity and the third component of the isospin {% = 0) are good quantum numbers, and the 
subscript a denotes collectively other non- fixed 11 quantum numbers. We use the j-j scheme 
for the coupling of the angular momenta, as illustrated in Fig. The Jacobi coordinates 
used in our study are depicted in Fig. [TJ This choice of coordinates allows us to separate 
the center-of-mass motion and guarantees that kinetic energy operator is independent of the 
angular variables. 

By projecting each of the Eqs. ([H]) on its natural configuration space basis one obtains 
coupled integro-differential equations for the radial FYAs yi,Zi)). Note that, contrary 

to the 3A^ problems, the number of radial FYAs is infinite even when the pair interaction is 
restricted to a finite number of partial waves. This divergence comes from the existence of 
the additional degree of freedom in the expansion of the K-tjpe components. Therefore 
an additional truncation is needed in numerical calculation by taking into account only the 
relevant amplitudes. 

C. Boundary conditions 

Eqs. (E]) needs to be supplemented with appropriate boundary conditions, which can be 
written in the Dirichlet form. For both bound and scattering states, the FYAs satisfy the 
regularity conditions: 

(0, y,, z,) = (x„ 0, z,) = (x„ y,, 0) = 0. (11) 

For bound state problems, since the wave functions are compact, the regularity conditions 
can be completed by forcing the amplitudes JF" to vanish at the borders of the hypercube, 

[0,^max] X [0, ^max] ^ [0,.^max]; 

i (^max) Vii ^i) j ("^i) ^^ax; ^i) ; (-^jj Vii ^max) 0' ('^^) 

The hypercube is chosen large enough to accommodate the wave functions. 

On the other hand, a scattering state near threshold contains two coupled channels, n-^He 
and p-'^H, both of which are of type K. For these amplitudes, we impose following matching 
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condition at Zi = ^max, 



2 y 5',vi^T-,i.«.T|e *"'iMj(r/,p;Z^ax)C7^3r,.,ii ) • (13) 



3 '2 2 



In the above equation, p„ = jmvrei is neutron's momentum of the n-^He channel, while Pp = 
[pn + 1^ (-SsH — -B^Hc)] is the proton's momentum of the p-^H recoils; hf"^ are the spherical 
Hankel functions, and ut is the outgoing Coulomb function of p-^B. channel with the strength 
7] = ^amc/ (hpp). The functions {f"''{xi,yi)} j^^^i j^^^j,^ with T| = ±| are the normalized 
Faddeev amplitudes for ^He and ^H, which we obtain by solving the corresponding 3A^ bound 
state problems. 

The wave functions are dominated by the total isospin T=0 state and, to describe the 
^He ground state, we can safely limit ourselves to the T=0 state. However, in order to 
separate the n-^He and p-^B. channels in the initial scattering state, we need to consider the 
admixture of total isospin T=l and T=2 states; we therefore solve the scattering problem 
without constraining the total isospin value. 

III. RESULTS 

In this work we have performed rigorous FY calculations for five sets of nuclear Hamil- 
tonians: Avl8, I-N3L0^N0Y, Avl8+UIX and I-N3L0+UIX*. Here Avl8 stands for the 
Argonne vlS potential I-N3L0 the chiral N^LO potential of Idaho group 0, INOY the 



semi-realistic configuration space potential that has been recently derived by Doleschall 20| , 



and UIX is the tri-nucleon interaction derived at Urbana [2ll]. For the case of I-N3L0+UIX*, 
we have slightly change the parameters in UIX so as to reproduce exactly the triton binding 
energy, see Ref. [ISi] for details. 

Our theory predictions for each of the model Hamiltonians used are listed in Table [B 
In addition to the hen cross section, the table also shows the binding energies of ^H, ^He 
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TABLE I: Theory predictions for the hen cross section, for five sets of the reahstic nuclear Hamil- 
tonians . We also list the values of binding energies (in units of MeV) for '^H, ^He and ^He, the 
point-proton rms radius ^] rHe4 (in fm) and the D-state probability PD(^He) (in per cents) of 
the a-particle as well as the spin-triplet n-^He scattering length anUeS (in units of fm) . The uncer- 
tainties attached to a represent the variation of a as the cutoff parameter A is varied in the range 
A = (500 ~ 900) MeV. 





be(3h) 


BE(3He) BE(''He) 




PD('He) 


OnHeS 


a [fib] 


Avl8 


7.623 


6.925 


24.23 


1.516 


13.8 


3.43 - 


0.0082i 


80.0 ± 12.2 


I-N3L0 


7.852 


7.159 


25.36 


1.52 


9.30 


3.56 - 


0.0070i 


57.3 ± 7.9 


INOY 


8.483 


7.720 


29.08 


1.377 


5.95 


3.26 - 


0.0058i 


34.4 ±4.5 


Avl8+UIX 


8.483 


7.753 


28.47 


1.431 


16.0 


3.23 - 


0.0054i 


49.4 ±8.5 


I-N3L0+UIX* 


8.482 


7.737 


28.12 


1.475 


10.9 


3.44 - 


0.0055i 


44.4 ±6.7 


Exp.: 


8.482 


7.718 


28.30 


1.475(6) 




3.278(53) 


-0.001(2)i 


55 ± 3, 54 ± 6 



and ^He, the point-proton rms radius rHe4 of "^He,*^ the D-state probability in "^He, and the 
spin-triplet n-^Re scattering length a„He3- 

The bound state properties calculated in this work agree well with those obtained in the 



other calculations 



24 



25 



26!]. Theoretical calculations for n-^He scattering are much less 
established, but we have checked that our results agree within 2% with the momentum-space 
FY calculation of the Lisboa group 27[ , as well as with the RGM calculation carried out by 
Hofmann 28| for the Avl8 and Avl8±UIX Hamiltonians. 

The table also indicates that the three-nucleon interactions (TNIs) play an important 
role in bringing the binding energies and the spin-triplet n-^He scattering length close to 
their respective experimental values. We remark that there is some uncertainty in the ex- 
perimental value of the scattering length. The value listed in the table is due to an i?-matrix 



analysis 



2^ based on the n-^He scattering data measured before the year 2002. Recently 



The point-proton rms radius rHo4 is defined as (rHo4)^ = ''c(^H*5) — — r'^, where rc('*He) is the proton 
charge rms radius of ^He, rp and r„ are the rms charge radius of the proton and neutron, respectively. 
See, for example, 22] for detailed explanation. With the ''He proton charge radius r-cCHe) — 1.681(4) fm 
obtained in a recent analysis 23| and with the 2008 PDG values for the proton and neutron rms radii, 
rp = 0.875(7) fm and rl = -0.1161(22) fm^ we arrived at met = 1.475(6) fm, which is about 1.4 % 
larger than the estimate given in [22] . 
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new measurements of the coherent scattering length have been performed at NIST SOj and 
ILL [31], but the results of the two groups do not agree with each other and both of them 
are in disagreement with the old ILL measurement. 

The D-state probability in the a-particle, P^i'^B.e), is closely related to the tensor forces 
and shows strong model dependence. However this quantity is not an observable, and it 
turns out to be difficult to constrain this quantity by studying other processes that are 
sensitive to Pol^He) [32 1. 

As can be seen from the Table, the calculated value of the hen cross section for Avl8 
(INOY) is too large (small), while the remaining three other Hamiltonians exhibit only a 
mild variation around the experimental values. The discrepancy of the Avl8 prediction is 
not surprising, since this Hamiltonian lacks TNIs and fails to reproduce even the binding 
energies. While a full discussion on model dependence will be given in the next subsection, 
this is a good place to discuss a feature that is specific to INOY. Although the INOY 
potential is capable of reproducing the binding energies and rms radii of A = 3 system 
quite accurately 33|], it is known to give too large binding energies and too small rms radii 
for A > 4 3 

and for nuclear matter [35|]. Such a feature leads to the reduction of the 
overlap integration between ra-^He and ^He and hence to the underestimation of the hen 
cross section. The INOY potential also deviates from the Tjon-line (a line that correlates 
the A = 3 and A = 4 binding energies) rather severely, which also indicates that caution 
should be exercised in using INOY for the A> A systems. 

For further discussion, we list in Table [Til the values of the hen matrix elements and the 
LECs evaluated at A = 700 MeV, for each of the nuclear Hamiltonians. The real part of Ai 
is written as the sum of IB and 2B contributions, with the dependence on g^g and (74^ also 
shown. Several remarks are in order here. Firstly, the 2B contributions turn out to be about 
two times as big as IB but with an opposite sign, which is consistent with the observation 



made in 



HI]. Secondly, there are substantial model-dependence even in the IB sector. 



which might be traced to the fact that not all the adopted Hamiltonians do not accurately 
reproduce the ERPs that govern the long-r contributions of IB. Finally, the inclusion of 
TNI(UIX) plays quite an important role in reducing the model dependence. 
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TABLE II: Theory predictions for M (in units of fm^/^), obtained with A = 700 MeV. is 
written in the format of: (IB) + (2B w/o CT) + {g^s-term) + (5f4^-term) = (total). 





9As 


94v 












Avl8 


0.3958 


0.1947 


0.1531 


-0.3777 + 0.0237 c/4s 


- 0.0403 g4v = 


-0.2231 


0.0249 


I-N3L0 


0.3919 


2.7479 


0.1304 


-0.2248 + 0.0198 54s 


- 0.0371 g4v = 


-0.1885 


0.0203 


INOY 


0.2313 


0.8021 


0.1045 


-0.2035 + 0.047354s 


- 0.0725 g4v = 


-0.1462 


0.0154 


Avl8+UIX 


0.5810 


-0.4615 


0.1518 


- 0.3567 + 0.0205 54s 


- 0.0377541, = 


-0.1756 


0.0179 


I-N3L0+UIX* 


0.5402 


2.3249 


0.1305 


-0.2253 + 0.017554s 


- 0.0347 541, = 


-0.1661 


0.0183 



A. Model-dependence 

In understanding the observed model dependence, it seems informative to recall the 



results of our previous study 13| of the Ml properties of the A < 3 nuclei. It was found 



m 



ISj that the Ml matrix elements (MEs) of the A=3 systems are strongly correlated with 
the triton binding energy B^. Since B3 governs the long-distance contributions to the MEs, 
the model-dependence (viz., variations in the MEs corresponding to the different nuclear 
Hamiltonians that give different values for ^3) cannot be cured by the renormalization 
procedure through the local (or short-ranged) operators. However, using the empirically 
established correlation curve between MEs and B3, we were able to calibrate the values of 
MEs by replacing the model-predicted B^ to its experimental value, which resulted in an 
essentially model-independent theoretical prediction. 

We expect that we can in principle adopt the same procedure for the hen process. Indeed 
we found that a is almost proportional to the quantity 

C=[q (flnHes/ rne4:f] ■ (14) 

We plot in Fig. [2] a versus C for the five nuclear Hamitonians under consideration. Fig. [2] 
suggests a strong correlation between a and If we take this correlation seriously and 



consider the quantity a defined by 



a ^ ^ a, (15) 



then a turns out to be almost model-independent: 

a = (54.8 ± 8.4, 54.5 ± 7.5, 53.9 ± 7.1, 54.8 ± 9.4, 56.4 ± 8.5) fib (16) 
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FIG. 2: The hen cross section a (in the unit of fib) plotted against C/Cexpj where ( = 



different choices of the cutoff parameter: A = 500 MeV (bottom), 700 MeV (middle with a filled 
circle) and 900 MeV (top). The lines in red with the label "Exp" denote the experimental data, 
54 ±6 Is'] and 55±3 jj]. 

for Avl8, I-N3L0, INOY, Avl8+UIX and I-N3L0+UIX*, all of which are in perfect agree- 
ment with the data. 

It is however not quite clear whether the correlation between a and ( is accidental or 
physical in nature. Furthermore, since hen is a four-body process that may involve a large 
number of ERPs, our numerical results that only cover five sets of Hamiltonians may not 
be sufficient to establish the correlation unambiguously. For example, a similar correlation 
behavior can be seen if we plot a against q^P'^^ . (It is however not clear whether this 
is a separate correlation or just another way to express the correlation between a and 
^, for Pd and rHe4 may not be independent of each other.) Without delving into the 
discussion of physics behind the correlation between a and we take here the viewpoint 
that the a calculated with a nuclear Hamiltonian that does not reproduce the relevant EPRs 
should be considered much less reliable than the a obtained with a nuclear Hamiltonian 
that does reproduce the relevant EPRs. Based on this viewpoint, we adopt here the results 



\_<1 (anHe3/'''He4)^] • For each Hamiltonian, the error bar denotes a range of variation for three 
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TABLE III: Cutoff-dependence of our prediction with the Avl8+UIX wave function. The values 
of g^s and g^y are determined so as to reproduce the magnetic moments of the '^H and '^He nuclei. 
For its dependence on g^g and g^^ are explicitly shown. Ai and a are in units of fm^/^ and 

fj,b, respectively. 



A [MeV] 


943 


94v 








a 


500 


0.8366 


1.9068 


-0.0915 + 0.0246^45 


-0.047154^ = 


-0.1608 


0.0166 


40.9 


600 


0.6990 


0.5886 


-0.1593 + 0.0231^45 


- 0.0433 c/4^ = 


-0.1686 


0.0173 


44.9 


700 


0.5810 


-0.4615 


-0.2049 + 0.0205^45 


- 0.0377 c/4^ = 


-0.1756 


0.0179 


48.7 


800 


0.4517 


-1.3622 


-0.2346 + 0.0176^45 


- 0.0319 g4^ = 


-0.1832 


0.0186 


52.9 


900 


0.3169 


-2.2069 


-0.2547 + 0.0149^45 


- 0.0265 g4v = 


-0.1915 


0.0195 


57.9 



obtained with Avl8+UIX and I-N3L0+UIX* as the most reliable predictions. For these 
two Hamiltonians, the calculated hen cross sections are: 

a = 49.4±8.5;u6 (Avl8 + UIX), a = 44.4 ± 6.7 /i6 (I-N3L0 + UIX*). (17) 

Both are lower by a few fib than the central value of the the experimental data but still in 
good agreement with the data within the error bars. 

B. CutofF dependence 

We now turn our attention to the cutoff dependence. Table IIIII shows to what extent 
our results for the Avl8+UIX wave functions change when the cutoff parameter A is varied 
over a range A = 500 ~ 900 MeV. Table IIIII indicates that the renormalization procedure of 
LECs {g4s and g^y) plays an essential role in removing the cutoff-dependence. As a way of 
quantifying its effectiveness, we define the quantity 

total (As) - A^total(Al) 
finite (As) - A^finitc(Al)' 

where the subscript "finite" stands for "finite-range term contributions" , and it corresponds 
to a case where all the terms other than the contact term (CT) contributions are included. 
Perfect renormalization invariance would correspond to -R=0. It turns out that Rhen = 
0.189 for (Ai, As) = (500, 900) MeV. Thus the renormalization procedure of LECs has 
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TABLE IV: Matrix elements for the Avl8+UIX wave function with A = 700 MeV; the LEG 



corresponding to this case 


are: g4,s = 0.581 and 


= -0.4615. 








H) 


/i( Hej 


)Kjvl 




IB: LO 


2.5727 


-1.7632 


0.0964 


-0.0136 


IB: RC 


-0.0171 


0.0037 


0.0554 


-0.0075 


IB-total 


2.5556 


-1.7595 


0.1518 


-0.0211 


2B: Ivr (NLO) 


0.2292 


-0.2258 


-0.1657 


0.0195 


2B: IvrC (N^LO) 


0.1578 


-0.1289 


-0.1465 


0.0172 


2B: 27r (N^LO) 


0.0419 


-0.0408 


-0.0445 


0.0052 


iilillc ILULcLi W / (J \^ L j 




— Z. lOOU 




W.UZUo 




0.0193 


0.0190 c/4s 


0.0205 gAs 


-0.0014(74. 




+0.0363 541, 


-0.0354 54^ 


-0.0377 


+0.0044 54„ 


2B: CT (N^LO) 


= -0.0055 


= 0.0274 


= 0.0293 


= -0.0029 


Total 


2.9790 


-2.1276 


-0.1756 


0.0179 



removed major part of cutoff-dependence, reducing the cutoff-dependence to about one-fifth 
of TVI finite- If is interesting to compare the above value of Rhen with the corresponding 
quantity Rhep obtained in a hep calculation in [sl. The hep calculation in [3J is based on 
the same MEEFT strategy (but used a different method for obtaining exact solutions to the 
nuclear Schroedinger equations). Ref. [3j reports Rhep = 0.137 for a slightly smaller window 
(Ai, A2) = (500, 800) MeV. Thus the previous hep calculation [3] is consistent with the 
hen calculation in the present work, and this consistency renders further support to the hep 
results in 



C. Details of ^A and chiral convergence 

Table HVl shows the individual contributions of the various IB and 2B terms to ;u(^H), 
/^(^He) and M/^en, calculated at A = 700 MeV for the Avl8+UIX Hamiltonian . We can 
see that the IB contribution to hen is highly suppressed due to the aforementioned orthog- 
onality between the initial and final wave functions. The NLO contribution, which comes 
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from the soft one-pion-exchange, is also suppressed for the Ml channel, due to the^ acciden- 
tal cancelation between the pion-pole and pion-seaguU diagram contributions [13]. These 
suppression mechanisms make chiral convergence rather unclear. For example, one might 
worry about the fact that the IttC and 27r contributions, both of which are N^LO, turn 
out to be comparable in size to the NLO Ivr contribution. It should be noted, however, 
that most of the IttC and 27r contributions are to be absorbed in the renormalization of the 
LECs, leaving very small effects on the net amplitudes. To demonstrate this, we define the 
effective matrix element, (^0^ effective^ ^ given operator O by 

^^^effective ^ ^^^^^^ _ ^^^^^^^ ^j^^^^^ ^^^^ _ (^g) 

Thus (^oy^'^'^^"'^ represents a change in the net amplitude that would occur if we omit the 
operator O. In evaluating the parenthesized quantity, the LECs should be readjusted so as 
to reproduce the experimental values of the A = 3 magnetic moments without (O); because 
of this readjustment we should expect ^(^^eSective _^ ^q^^ 

(^74s,^?4.) = (0.5223, 5.8848), 
igis,9iv) = (1.2433, 3.53455), 
(^74s,^?4.) = (0.5843, 0.6921), (20) 

where we have also listed the corresponding values of the LECs, which should be compared 
with ((74s, (yf4t,) = (0.581, —0.4615) that corresponds to the full calculation up to N'^LO. Eq. fl20|) 
demonstrates that the effective contributions of IvrC and 27r are very small, only about 6 
% and 2 %, respectively, relative to the values one would naively expect. This is in sharp 
contrast with the NLO soft one-pion-exchange, whose contribution cannot be absorbed in 
the LECs. A rigorous proof of chiral convergence would require a calculation that goes one 
order higher than considered in the present work; i.e., we need to go up to N^LO. 

IV. DISCUSSION AND CONCLUSION 

There can also be a "completely" consistent approach where both nuclear Hamiltoni- 
ans and transition operators are obtained in the same EFT framework. This approach 
requires much more involved calculations than MEEFT, but recent significant progress in 
constructing EFT Hamiltonians makes it more attractive. Also available is a pionless EFT 
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^^^^effective ^ 0.0749 - 0.0087z, 

^^^^^effective _ -0.0093 + 0.0004z, 

^2^^cffcctive _ -0.0010 + 0.000n. 



approach [37j where the matrix elements are evaluated perturbatively, but, unless it is ca- 
pable of reproducing all the relevant ERPs of the nuclear systems under consideration, its 
usefulness is limited. 

In this work we have performed an ab initio parameter-free calculation for the hen cross 
section, with the EM currents that have been derived from HBChPT up to N^LO. Exact 
nuclear wave functions have been obtained by solving the Faddeev-Yakubovski equations for 
realistic nuclear Hamiltonians. The calculated value of a shows a high degree of stability as 
the cutoff parameter A is varied over a wide range, A = (500 ~ 900) MeV, and we obtain 
as teh best estimate a = 49 A ± 8.5 fib for Avl8+UIX and 44.4 ± 6.7 fib for I-N3L0+UIX*. 
These values are in good agreement with the data, 54 ± 6 [5] and 55 ± 3 6|. 

The successful application of MEEFT to hen renders strong support to the previous 
MEEFT calculation for hep in Ref. [s]; furthermore, it demonstrates the basic soundness 
of the MEEFT strategy in general. The present treatment is open to several improvements 
such as: the inclusion of the next order in chiral perturbation, in particular the incorporation 
of the three-nucleon currents; a more stringent control of mismatch in the chiral counting 
between SNPA and a formally accurate chiral expansion that enters in the currents; a better 
understanding of the role that the counter terms play in the renormalization group property, 
etc. It is reasonable, however, to expect that the effects of these improvements are essentially 
accommodated in the above-quoted error estimate based on the cutoff dependence. A robust 
estimation of the hep S-factor has been a long-standing challenging problem in nuclear 
physics 3^. We believe that our MEEFT calculations of hep and hen have solved this 
problem to a satisfactory degree. Short of doing a fully formally consistent EFT calculation, 
which is at present out of our reach, MEEFT seems to be the best currently available 
approach. 
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